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ABSTRACT 

In the double detonation scenario for Type la supernovae it is suggested that a detonation initiates 
in a shell of helium-rich material accreted from a companion star by a sub-Chandrasekhar-mass White 
Dwarf. This shell detonation drives a shock front into the carbon-oxygen White Dwarf that triggers a 
secondary detonation in the core. The core detonation results in a complete disruption of the White 
Dwarf. Earlier studies concluded that this scenario has difficulties in accounting for the observed prop- 
erties of Type la supernovae since the explosion ejecta are surrounded by the products of explosive 
helium burning in the shell. Recently, however, it was propo sed that detonations might be possible 
for much less massive helium shells than previously assumed ([Bildsten et al.ll2007l) . Moreover, it was 
shown that even det onations of these minimum helium shell masses robustly trigger detonations of the 
carbon-oxygen core (|Fink et al.ll2"0 10). Therefore it is possible that the impact of the helium layer on 
observables is less than previously thought. Here, we present time-dependent multi-wavelength radia- 
tive transfer calculations for models with minimum helium shell mass and derive synthetic observables 
for both the optical and 7-ray spectral regions. These differ strongly from those found in earlier sim- 
ulations of sub-Chandrasekhar-mass explosions in which more massive helium shells were considered. 
Our models predict light curves which cover both the range of brightnesses and the rise and decline 
times of observed Type la supernovae. However, their colours and spectra do not match the obser- 
vations. In particular, their B — V colours are generally too red. We show that this discrep a ncy is 
mainly due to the composition of the burning products of the helium shell of the iFink et all ([20101 ) 
models which contain significant amounts of titanium and chromium. Using a toy model, we also 
show that the burning products of the helium shell depend crucially on its initial composition. This 
leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed 
Type la supernovae may still be feasible but further study of the shell properties is required. 
Subject headings: methods: numerical - radiative transfer - supernovae: general 



1. INTRODUCTION 

Although the standard single-degenerate 

(SD) Chandrasekh a r-mas s scenario (see 

iHillebrandt fc Niemeverl 120001 for a review) is capa- 
ble of explaining most of the o bserved diversity of 
Type la supernov ae _(SNe la) ([Hoflich fe Khokhlovl 
199 61 iKasen et all [2009) via the delayed-detonation 
(|KhokhlovllT99ll )" model, it suffers from severe problems 
in explaining the observed rate of SNe la. In particular, 
binary evolution population synthesis calculations 
predict rates which are an order of magni tude too low 
comp ared to the observed rate of SNe la ([Ruiter et al.1 
I2009L but see iMeng fc Yandl2010l ). In addition, recent 
observational studies suggest that Chandrasekhar-mass 
explosions of hydrogen-accreting carbon-oxygen (C/O) 
White Dwarfs (WDs) in SD binary systems cannot 
account for all SNe la: iGilfanov fc Bogdanl (pOlOl ) 
found that the X-ray flux of nearby elliptical galaxies is 
significantly weaker than expected for a population of 
WDs accreting hydrogen towards the Chandrasekhar- 
mass needed to explain th e observed super nova rate in 
elliptical galaxies (see also[Di Stcfano 2010). Moreover, 
there is growing observational e vidence that there are 
different populations of SNe la (jMannucci et al.1 120051 : 
iScannapieco fc Bildstenl[2005[ ). 

This has led to a revived interest in alternative explo- 
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sion mechanisms. Here we consider the double detona- 
tion scenario applied to sub-Chandrasekhar-mass C/O 
WDs. In that scenario a helium-accreting C/O WD 
explodes below Mch due to a detonation in the ac- 
creted helium shell which triggers a secondary core det- 
onati o n by compressional h e ating dWooslev fc Weaverl 
Il99l ILivne fc ArnetH [l995t IFink et all 12007^ This 
model has some very appealing features. Depending on 
the initial mass of the WD, a wide range of explosion 
strengths can be real i zed fe.g. [Wooslev fc Weaverl 1 1994t 
ILivne fc Arnettl 1 19951: IHoflich fc Khokhlovl 1 19961). More- 
over, population synthesis studies ([Ruiter et al.l I2009D 
predict rates comparable to the observed galactic super- 
nova rate. 

However, earlier work (IHoflich fc Khokhlovl 119961 : 
IHoflich et al.ll996HNugent et al.lll997t) found light curves 
and spectra of such models to be in conflict with the ob- 
served spectra and light curves of SNe la. The differences 
were mainly attributed to the composition of the outer 
layers. Due to the initial helium detonation in the outer 
shell, the ejecta of sub-Chandrasekhar-mass models are 
surrounded by a layer of helium and its burning products 
(which can include iron-peak nuclei). This, however, is in 
apparent contradiction to the layered composition struc- 
ture of observed SNe la, where the composition changes 
from iron-group elements in the core to lower mass ele- 
ments in the outer layers. 

In a preceding paper iSim et all ([2010) have shown 
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that artificial explosions of "naked" sub-Chandrasekhar- 
mass WDs can reproduce the observed diversity of 
SNe la. Thus it is natural to ask if somewhat mod- 
ified properties in the initial helium shells of realistic 
sub-Chandrasekhar-mass models can reduce the nega- 
tive post-expl osion effect of th i s shel l on the observables. 
In particular, iBildsten et al.l (|2007l ) recently presented 
new calculations, indicating that detonations might oc- 
cur for m uch less ma s sive h elium shells than previously 
thought. iFink et al.l (12010ft adopte d the minimum he- 
lium shell masses of IBildsten et al.l (|2007l ) and investi- 
gated whether such low-mass helium detonations are ca- 
pable of triggering a secondary detonation in the C/O 
core of the WD. In that study, they concluded that 
as soon as a detonation in the helium shell initiates, a 
subsequent core detonation is virtually inevitable. For 
example, they found that even a helium shell mass as 
low as 0.039 Mq is sufficient to detonate a C/O WD of 
1.125 M Q . 

Here, we focus on the obser v able p roperties of the 
models presented in IFink et ail (|20100 and their com- 
parison to real SNe la. In particular, we investigate 
whether the low helium shell masses of these models help 
to alleviate the problems encountered previously when 
comparing double detonation sub-Chandrasekhar-mass 
models to observed spectra and light curves of SNe la 
(Hofl ich fc Khokhlovlll99l INueent et allll99l . 

The outline of the paper is as follows: i n Section [2] 
we giv e a short summary of the models of IFink et al.1 
(2010) before briefly describing details of our radiative 
transfer simulations in Section al In Secti o n [1] we present 
synthetic observables for the IFink et all (|2010D models 
and compare them to the observed properties of SNe la. 
The results of this comparison and implications for future 
work on sub-Chandrasekhar-mass double detonations are 
discussed in Section [5] Finally, we draw conclusions in 
Section [5] 

2. MODELS 

Adopting the minimum helium shell masses required 
to initiate a helium detonation in the s hell according to 
IBildsten et all (|2007f ). IFink et al.l (|2010f ) investigated the 
double detonation scenario for six models representing a 
range of different C/O core masses (the model param- 
eters are summarized in Table [1]). In all their models, 
they ignite an initial helium detonation in a single point 
at the base of the helium shell located on the positive 
z-axis (hereafter referred to as the "north-pole" of the 
WD). From there the helium detonation wave sweeps 
around the WD core until it converges at the south pole. 
At the same time a shock wave propagates into the core 
and converges at a point off-centre. Finding conditions 
that might be sufficient to initiat e a core detonatio n in a 
finite volume around this point, IFink et ail (|2010| ) then 
trigger a secondary core detonation at that point. This 
secondary detonation disrupts the entire WD and yields 
ejecta with a characteris tic abundance dist ribution which 
is shown for Model 3 of IFink et ail (pOlOl ) in Figure [TJ 

In the initial helium shell the burning does not reach 
nuclear statistical equilibrium due to the low densities. 
Thus it mainly produces iron-group elements lighter than 
56 Ni such as titanium, chromium and iron - includ- 
ing some amount of the radioactive isotopes 48 Cr and 
52 Fe. Aside from a small mass of calcium, no significant 



amounts of intermediate-mass elements are produced in 
the shell. However, a large fraction of helium remains 
unburned. 

The core detonation yields both iron-group and 
intermediate-mass elements, the relative amounts of 
which depend crucially on the core density and thus 
WD mass. Models with more massive WDs produce 
more iron-group material and less intermediate-mass el- 
ements (similar to the explosions of the naked sub - 
Chandrasekhar-mass WDs studied in iSim et alH 2010). 
The m ost massive Model 6 of the sequence of lFink et al.l 
(2010) produces almost only iron-group material and 
hardly any intermediate-mass elements (the most impor- 
tant burning products in the shell and core of the differ- 
ent models are listed in Table [1]). 

Due to the single-point ignition the models show strong 
ejecta asymmetries which can be divided into two main 
categories. First, the helium shell contains more iron- 
gr oup materia l on th e northern hemisphere as discussed 
by IFink et al.l (|2010D . Second, the ignition point of the 
secondary core detonation is offset from the centre-of- 
mass of the model due to the off-centre convergence of 
the shock waves from the helium detonation. 

3. RADIATIVE TRANSFER SIMULATIONS 

To derive synthetic observables for the models, we 
performed radiative transfer simulations with the time- 
dependent multi-dimensional Monte Carlo radiative 
transfer code ARTIS described by iKromer fc Siml (2009) 
and Urn] (pfffll . Since the models produce significant 
amounts of 48 Cr and 52 Fe, we extended ARTIS to take 
into account the energy released by the decay sequences 
52 Fe -> 52 Mn -» 52 Cr and 48 Cr -» 48 V -» 48 Ti in ad- 
dition to the radioactive decays of 56 Ni — > 56 Co and 
56 Co — » ■ 56 Fe which form the primary energy sou rce of 
SNe la (|Truran et al.lll967L ICokate fc McKeelil969fl . For 
other radioactive nuclei which are synthesized during the 
explosion in a non-negligible amount (e.g. 44 Ti in the 
shell) the life times are much longer than those of the 
56 Ni decay-sequence. Thus they can be neglected at early 
times when the decays of 56 Ni and 56 Co power the light 
curves. 

The total 7-ray energy E to t that will be emitted from 
t = to t — > 00 in the decay chains we consider, will be 

-Etot = (Ese m + E'oBQo)— — + (E52 Fc + E52 Mn )— — 

777,56 N; 77152 Fe 

+ (£J tICr + .E«v)— (!) 

777,480 

Here Esem, (Es6q , i?52p c , i?52M n , E4sq t , i?48y) is the 
mean energy emitted per decay of 56 Ni, ( 56 Co, 52 Fe, 
52 Mn, 48 Cr, 48 V). Similarly, M5e Ni (M52 Fe , M4s Cr ) is 
the initial mass of 56 Ni ( 52 Fe, 48 Cr) synthesized in the 
explosion and 77756 Ni (m52 Fc , m4s Cr ) the mass of the 56 Ni 
( 52 Fe, 48 Cr) atom. 

Following EHci J2005), this energy is quantized into 
J\f = E to t Ao identical indivisible energy "pellets" of ini- 
tial co- moving frame (cmf) energy eo. The pellets are 
first assigned to one of the decay sequences in propor- 
tion to the amount of energy deposited in the different 
decay sequences (terms on the right-hand-side of Equa- 
tion [lj. Then they are distributed on the grid accord- 



Synthetic observables for sub-Chandrasekhar-mass models of SNe la 
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Table 1 

Parameters of our model sequence. 



Model 


1 




2 




3 




4 




5 




6 




3c 




3 in 




Mtot 


936 




1.004 




1.080 




1.164 




1.293 




1.3885 




1.025 




1.071 




-Aicorc 


0.810 




0.920 




1.025 




1.125 




1.280 




1.3850 




1.025 




1.025 




M core ( 56 Ni) a 


1.7 x 10" 


-l 


3.4 x 10" 


i 


5.5 x 10" 


l 


7.8 x 10" 


l 


1.05 




1.10 




5.5 x 10" 


-l 


5.6 x 10" 


l 


M core ( 52 Fe) a 


7.6 x 10" 


-3 


9.9 x 10" 


-3 


9.6 x 10" 


3 


7.9 x 10" 


3 


4.2 x 10" 


3 


1.7 x 10" 


4 


9.6 x 10" 


-3 


9.4 x 10" 


3 


M core ( 48 Cr) a 


3.9 x 10" 


-4 


4.6 x 10" 


4 


4.5 x 10" 


4 


3.8 x 10" 


4 


2.1 x 10" 


4 


7.1 x 10" 


5 


4.5 x 10" 


4 


4.4 x 10" 


4 


M rorr ( 44 Ti) 


7.2 x 10" 


-6 


9.0 x 10" 


(i 


1.1 x 10" 


5 


1.4 x 10" 


5 


1.4 x 10" 


5 


9.9 x 10" 


6 


1.1 x 10" 


-5 


1.2 x 10" 


5 


M rorf ,( 40 Ca) 


2.0 x 10" 


-2 


2.1 x 10" 


2 


1.8 x 10" 


2 


1.4 x 10" 


2 


6.9 x 10" 


3 


1.8 x 10" 


3 


1.8 x 10" 


-2 


1.8 x 10" 


2 


Mcorc( 36 Ar) 


2.2 x 10" 


-2 


2.2 x 10" 


-2 


1.9 x 10" 


2 


1.5 x 10" 


2 


6.7 x 10" 


3 


1.7 x 10" 


3 


1.9 x 10" 


-2 


1.9 x 10" 


2 


M core ( 32 S) 


1.3 x 10" 


-1 


1.2 x 10" 


1 


1.0 x 10" 


1 


7.5 x 10" 


2 


3.3 x 10" 


2 


8.0 x 10" 


3 


1.0 x 10" 


-1 


1.0 x 10" 


1 


M ro rrf 28 Si) 


2.7 x 10" 


-1 


2.5 x 10" 


1 


2.1 x 10" 


1 


1.4 x 10" 


1 


6.1 x 10" 


2 


1.5 x 10" 


2 


2.1 x 10" 


-1 


2.0 x 10" 


1 


Mcore( 24 Mg) 


4.5 x 10" 


-2 


3.5 x 10" 


-2 


2.4 x 10" 


2 


1.1 x 10" 


2 


9.3 x 10" 


3 


4.3 x 10" 


3 


2.4 x 10" 


-2 


2.2 x 10" 


2 


M C ore( 16 0) 


1.4 x 10" 


-1 


1.1 x 10" 


1 


8.0 x 10" 


2 


4.2 x 10" 


2 


3.1 x 10" 


2 


1.2 x 10" 


2 


8.0 x 10" 


-2 


7.7 x 10" 


2 


M coro ( 12 C) 


6.6 x 10" 


-3 


4.4 x 10" 


-3 


2.7 x 10" 


3 


8.8 x 10" 


4 


5.9 x 10" 


3 


7.4 x 10" 


4 


2.7 x 10" 


-3 


3.3 x 10" 


3 


Af sh 


0.126 




0.084 




0.055 




0.039 




0.013 




0.0035 








0.046 




M sh ( 56 Ni) a 


8.4 x 10" 


-4 


1.1 x 10" 


-3 


1.7 x 10" 


3 


4.4 x 10" 


3 


1.5 x 10" 


3 


5.7 x 10" 


4 


- 




1.1 x 10" 


8 


Af, h ( 52 Fe) a 


7.6 x 10" 


-3 


7.0 x 10" 


-3 


6.2 x 10" 


3 


3.5 x 10" 


3 


1.2 x 10" 


3 


2.0 x 10" 


4 






6.1 x 10" 


9 


AL h ( 48 Cr) a 


1.1 x 10" 


-2 


7.8 x 10" 


-3 


4.4 x 10" 


3 


2.2 x 10" 


3 


6.8 x 10" 


4 


1.5 x 10" 


4 






7.9 x 10" 


8 


/If t. C 44 Ti"l 
iw sh1, rij 


7.9 x 10" 


-3 


5.4 x 10" 


-3 


3.4 x 10" 


3 


1.8 x 10" 


3 


4.9 x 10" 


4 


6.2 x 10" 


5 






4.2 x 10" 


5 


Jw shl ^ a J 


A 7 v 1 n~ 


-3 


q o v in - 

G.Z X 1U 


-3 


9 n v 1(1- 


3 


Q 9 y 1(1- 


3 


O.o X 1U 


4 


9 i v 1 


4 






9 CI y 1(1- 

O.O X 1U 


3 


JW s h( Arj 


5.0 x 10" 


-3 


3.2 X 10" 


-3 


2.0 x 10" 


3 


1.6 X 10" 


3 


5.2 x 10" 


4 


1.3 X 10" 


4 






2.7 x 10" 


2 




2.2 x 10" 


-3 


1.2 x 10" 


-3 


7.8 x 10" 


4 


1.3 x 10" 


3 


4.3 x 10" 


4 


1.9 x 10" 


4 






1.4 x 10" 


2 


M sh ( 28 Si) 


4.8 x 10" 


-4 


2.5 x 10" 


4 


1.4 x 10" 


4 


4.7 x 10" 


4 


1.6 x 10" 


4 


1.3 x 10" 


4 


— 




5.9 x 10" 


4 


M sh ( 24 Mg) 


4.3 x 10" 


" 


2.3 x 10" 


-5 


1.3 x 10" 


5 


1.0 x 10" 


4 


3.8 x 10" 


5 


3.6 x 10" 




— 




3.0 x 10" 


5 


M sh ( 16 0) 


3.2 x 10" 


-6 


1.7 x 10" 


(i 


1.9 x 10" 


6 


6.0 x 10" 


5 


2.1 x 10" 


5 


2.6 x 10" 


5 


- 




1.0 x 10" 


5 


M sh ( 12 C) 


1.2 x 10" 


-3 


5.3 x 10" 




2.2 x 10" 




7.9 x 10" 


5 


1.7 x 10" 


5 


1.6 x 10" 


(i 






3.2 x 10" 


5 


Af sh ( 4 He) 


8.3 x 10" 


-2 


5.3 x 10" 


-2 


3.3 x 10" 


2 


2.0 x 10" 


2 


6.9 x 10" 


3 


1.7 x 10" 


3 






2.7 x 10" 


4 


Ami5(i() / mag 


0.88 




1.25 




1.74 




1.77 




1.23 




1.29 




1.63 




1.37 




tmax(B) / days 


17.0 




17.0 




17.7 




16.4 




15.2 




14.3 




18.2 




18.1 




AfB, max / mag 


-15.9 




-17.3 




-18.4 




-19.3 




-19.8 




-19.9 




-19.2 




-19.1 




My, max 1 mag 


-17.6 




-18.8 




-19.6 




-19.9 




-20.1 




-20.1 




-19.4 




-19.4 




AfR, max / mag 


-18.4 




-19.1 




-19.4 




-19.4 




-19.4 




-19.2 




-19.0 




-19.1 




Mi, max / mag 


-18.9 




-19.2 




-19.2 




-19.4 




-19.6 




-19.7 




-18.9 




-19.0 




(U - B) Bm ax / mag 


0.58 




0.41 




0.50 




0.19 




0.08 




-0.26 




-0.04 




0.35 




(B - V) B max / mag 


1.67 




1.51 




1.13 




0.59 




0.28 




0.08 




0.11 




0.30 




(V - i?)Bmax / mag 


0.74 




0.36 




-0.13 




-0.46 




-0.67 




-0.82 




-0.30 




-0.28 




(V - 7)Bmax / mag 


1.26 




0.39 




-0.47 




-1.05 




-1.48 




-1.47 




-0.64 




-0.49 





Note. — Aftot j Af core , and Af s h are the masses of the WD, the C/O core, and the helium shell, respectively. All masses arc given in 
units of the solar mass. Amis(B) and i max (-B) refer to the decline parameter and the rise time to maximum light in the anglc-avcragcd 
.B-band light curves, respectively. Mb, max, Afv,max, max and Afi. max denote the anglc-avcragcd peak magnitudes at the true peaks 
in the given bands. Colours are quoted at time [t m ax(S)] of B-band maximum. 

a At maximum light 56 Ni, 52 Fe and 48 Cr will have mostly decayed to 56 Co, 52 Cr and a mixture of 4S V and 48 Ti, respectively. 



og (Mass fraction) 



-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 




10 20 10 20 10 20 10 20 10 20 10 20 



v in 1 3 km s 



Figure 1. Final composition structure of selected species of Model 3 of Fink et al. (2010). The individual panels show the mass fractions 
of He, O, Si, Ca, Ti and Ni (from left to right, respectively). The model is radially symmetric about the z-axis. 
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ing to the initial distribution of 56 Ni, 52 Fe and Cr, as 
appropriate, and follow the homologous expansion until 
they decay. Decay times are sampled in a two step pro- 
cess. If the pellet was assigned to the 56 Ni chain, for 
example, we first choose whether it belongs to a decay 
of the parent nucleus 56 Ni or the daughter nucleus 56 Co 
by sampling the probabilities Esb^/ (Ese^i + E^q q ) and 
E56q / (-E56N; + EssQg), respectively. Finally an appro- 
priate decay time is sampled 

idecay ( 56 Ni) = T ( 56 Ni) log (z X ) (2) 
idecay ( 56 Co) = -T ( 56 Ni) log {z X ) - T ( 56 Co) log (z 2 ) 

(3) 

from the mean life times of the 56 Ni [r ( 56 Ni) = 8.80 d] 
and 56 Co [r ( 56 Co) = 113.7 d] nuclei, m are random 
numbers between and 1. The 52 Fe [r ( 52 Fe) = 0.4974 d, 
r( 52 Mn) = 0.02114 d] and 48 Cr [r ( 48 Cr) = 1.296 d, 

r ( 48 V) = 23.04 d] chains are treated in the same way. 

Upon decay, a pellet transforms to a single 7-packet 
representing a bundle of monochromatic-radiation of cmf 
energy eo- The cmf photon energy of the 7-packets 
is randomly sampled from the relative probabilities of 
the 7-lines in the appropriate decay of the selected 
decay sequence - including annihi l ation lines due to 
positron emission. Following iLucvl ()2005l ). we assume 
that positrons released by radioactive decays annihilate 
in situ, giving rise to the emission of two 511 keV 7- 
ray photons. In doing so, we neglect the kinetic energy 
released by stopping the positrons, any positron escape 
and possible positronium formation which gives rise to 
the emission of c ontinuum photons (see discussion by 
iMilne et aLll2004D . Thus, our prediction of the 511 keV 
line flux should be considered as an upper limit. The 
7-packets are then propag ated through the ejecta as de- 
scribed bv lKromer fc Siml (|2009f ). 

The 7-line data i s take n from Table 1 of 
lAmbwani fc Sutherland! (|1988l) for t he 56 Ni decay- 
sequence and from iBurrowa (|2006h for the 48 Cr 
decay-sequence, respectively. Owing to the compara- 
tively short life times of 52 Fe and 52 Mn, these nuclei 
have mostly already decayed to their daughter nuclei 
when we start the radiative transfer simulation at ~ 1 
day. Since the ejecta at these early times are almost 
opaque to 7-rays, the 7-packets released by the 52 Fe 
decay chain will be thermalized rapidly. Therefore we 
do not follow the propagation of the 7-packets released 
by the 52 Fe decay chain, but immediately convert their 
energy to thermal kinetic energy (fc-packets in the 
framework of ARTis). 

The input models (composition, density, velocities) 
for the radiative transfer simulations were derived us- 
ing the trace r particles from the hydrodynamics sim- 
ulations (see iFink et ail 120101 for details on the hydro 
setup and use of tracer particles for nucleosynthesis). 
Since the tracer particles are Lagrangian, we use an ap- 
proach similar to the reconstruction of the density field 
in smooth particle hydrodynamics (SPH) simulations to 
construct the input model from the tracers. For the mod- 
els described here, the densit y field was obtained using 
the SPH method described bv lDolag fc Stasvszvnl (2009) 
(specifically, their equations 1-3) adopting N — 32 for 
the SPH smoothing-length normalisation factor. For the 



centre (xi) of each grid cell (i) in the model, the mass 
fractions (Xz,i) of the elements considered (Z = 1 to 30) 
were reconstructed using 

X z ,i = pi^Wdxi - .i-j .h,)Xx.,. (4) 

3 

where pi is the reconstructed mass density, Xz,j is the 
mass-fraction of element Z for tracer particle j (which 
lies at position Xj), hj is the SPH particle smoothing 
length and W{x, h) the SPH kernel function (de fined via 
equations 3 and 1 of iDolag fc~S tasvszvn 200^, respec- 
tively). The mass fractions of the important radioactive 
isotopes ( 56 Ni, 56 Co, 52 Fe and 48 Cr) were reconstructed 
from the tracer particle yields in exactly the same man- 
ner. The reconstruction was performed on an 80 x 160 
(r, z)-grid using the final state of the tracer particles at 
the end of the simulations (which were run up to the 
phase of homologous expansion) . 

For the radiative transfer simulation this 2D model is 
re-mapped to a 3D Cartesian grid of size 100 3 which 
co-expands with the ejecta. We then follow the radia- 
tive transfer from 2 to 120 days after the explosion, dis- 
cretized into 111 logarithmic time steps each of duration 
Aln(i) = 0.037. Using our detailed i onization treat- 
ment and the cd23_gf-5 atomic data set (|Kromer fc Sirnl 
2009), we simulated the propagation of 2 x 10 7 pack- 
ets. To speed up the initial phase of the calculation we 
made use of our initial grey approximation for the first 
30 time steps (adopting a parameterised grey opacity for 
the highly optically thick cells in the centre of our simu- 
lation volume). The first ten time steps were treated in 
local thermodynamic equilibrium (LTE). 

4. SYNTHETIC OBSERVABLES 

In this S ection we p r esent synthetic observables for the 
models of IFink et all (|2010D . First we consider angle- 
averaged ultraviolet, optical and infrared light curves and 
explore the diversity within this set of models and com- 
pare it to that of observed SNe la (Section 14. ip . Then 
we investigate the colour evolution and spectra of these 
models (Sections 14.21 and I4.3|) . In Section |4~41 we study 
the effects of the asymmetric ejecta composition. Finally, 
in Section 14.51 we briefly discuss the 7-ray emission from 
these models. 

4.1. Broad-band light curves 

In Figure [5] we show the angle-averaged ultraviolet- 
optical-infrared (UVOIR) bolometric and U, B. V, 
R, I, J, H, and K ba nd-limited li ght c urves for 
the model sequence of IFink et all (]2010f ) as ob- 
tained from our radiative transfer simulations. For 
comparison to observ ations we included ph otometric 
data of SNe 2005cf (iPastorello et all l2007bl) . 2 004eo 
(|Pastorello et al.ll2007al). 2001e l (iKrisciunas et al.ll2003T) 
and 1991bg (|Filippenko et all Il992t iLeibundgut et alj 
I1993D in the figure. While SN 1991bg, the prototypi- 
cal event of the subluminous 1991bg-like objects, marks 
the faint end of observed SNe la, SNe 2005cf, 2004eo and 
2001el are representative of the spectroscopically normal 
objects. 

Since our models form a sequence of increasing 56 Ni- 
mass, the peak brightness of the synthetic light curves 
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Figur e 2. Angle-averaged UVOIR bolometric and U \B,V ,R,I ,J \H ,K band-limited light curves for the model sequence of IFink et al.l 
(120101 ) as indicated by the colour coding. For compa rison photometr ical data of the spectroscop ically normal SN e 2005cf (triangles; 



Pastorcllo ct al. 2007b), 2004co (upside-down triangles; Pastorcllo ct al. 2007a) and 2001el (squares; Krisciunas ct al. 2003) as well as for 
the subluminous SN 1991bg (circles; Filippenko et al. 1992, Lcibundgut et al. 1993) are shown. 
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Figure 3. Angle-averaged colour curves of our models as indicated by the labels. The different panels correspon d to B — V, V - 
V — I colour (from left to right). Fo r comparison, the colou rs of the spectroscopically normal SNe 2005cf (triangles; Pastorcllo ct al. 



and 2004eo (upside-dow n triangles; [Pastorcllo ct al. 2007aj) as well as for the the subluminous SN 1991bg (circles; Filippenko ct a 
Leibundgut ct al. 1993) are shown. 
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increases from Model 1 to Model 6: for B-band max- 
imum, for example, we find a range from —15.9 mag 
for Model 1 to —19.9 mag for Model 6 (values for the 
intermediate models and other bands are given in Ta- 
ble [l}. This covers almost the full r ange of peak bright - 
nesses observed among SNe la (e.s. IHicken et~al~l l2009 N >. 
excluding only the very brightest events like SN 2003fg 
(|Howell et al.ll2006h which have been suggested to derive 
from su per-Chandrasekhar-mas s progenitor syste ms (see 
however IfflTlebrandt et aTll2007t iSim et al.ll2007h . 

Despite their lower total mass, and thus lower overall 
opacity, our model light curves do not show particularly 
fast evolutio n compared to Chandrasekh ar-mass models. 
In contrast, [Hoflich fc Khokhlovj (|1996| ) reported a fast 
initial rise and broad peaks for their helium detonation 
models. For the rise times from explosion to B-band 
maximum [t max (B)] we find values between 17.0 and 17.7 
days for our least massive models (1 to 3, cf. Tabled]). 
This is in good agreem ent with observational findings 
(e.g. IHavden et afll2010l find rise times to S-band max- 
imum between 13 and 23 days with an average value of 
17.38 ±0.17 days). In contrast, the more massive Models 
(4 to 6) show shorter i?-band rise times with increasing 
mass (see Tabic [TJ): for Model 6, t max (B) is only 14.3 
days. The peak time is mainly set by the diffusion time 
for photons to leak out from the 56 Ni-rich inner core. Ow- 
ing to the larger densities in the more massive models, 
nuclear burning produces 56 Ni out to much higher veloc- 
iti es than in the lo wer mass models (compare Figure 3 
of lFink et alJ l2010). Moreover, the helium shell masses 
in our models decrease for the heavier WDs. Thus, the 
mass (and therefore opacity) on top of the 56 Ni-rich in- 
ner core decreases with increasing WD mass and photons 
start to escape earlier for the more massive models mak- 
ing their light curves peak faster (compare also the rise 
times to bolometric maximum in iFink et al.ll2010h . In 
addition, we note that the radioactive nuclides which are 
produced in the helium shell burning have some influence 
o n the initial rise p hase of the light curves (see Figure 8 
of IFink et aLll2010h . 

In the post-maximum decline phase, characterized by 
the decline parameter Atoi5(_B) which gives the change 
in B-band magnitude between maximum and 15 days 
thereafter, our models show a peculiar characteristic. 
Observationally, brighter S Ne la show a trend of slower 
decli ning light curves (e.g. iPhillipsI [1991 IHicken et"all 
2009). In contrast, the model light curves decline faster 
along our sequence from Model 1 to 4 despite increasing 
brightness (cf. Table [I). 

Specifically, we find Ami 5 (B) ~ 0.88 for Model 1, 
which according to peak brightness would be classi- 
fied as a subluminous explosion. However, observa- 
tionally these events are chara cterized by a fast de- 
cline \Am^(B) ~ 1.9 e.g. iGarnavich et al.1 120041 : 
iTaubenberger et all 120081 ] . In contrast, Model 4, which 
has a B-band peak magnitude typical for a spectroscop- 
ically normal SN la, yields Ami^(B) ~ 1.77 which is 
much faster than t ypically observed [Amis(S) ~ 1, e.g. 
IHicken et al.l 12009} . Although Amig(B) is lower for the 
brighter Models 5 and 6, they still decline too fast com- 
pared to observed objects of corresponding brightness. A 
similar trend for more rapidly declining light curves from 
more massive models is also seen in t he bolometric ligh t 
curves of Model 2 to 6 (cf. Table 4 of lFink et af1 [20101. 



In contrast to earlier work by IHoflich fc Kho khlov 
(1996), who studied sub-Chandrasekhar-mass models 
with significantly more massive shells, we do not find 
particularly fast evolution in the post-maximum de- 
cline of our light curves compared to those obtained for 
other types of models. Applying our radiative transfer 
code ARTIS for example to the well-known W7 model 
(|Nomoto et al J 119841 : iThielemann et alJ ll986). which is 
widely regarded as a good standard for SNe la, we find 
Ami§(B) ~ 1.6 when using the same atomic data set as 
adopted here. This value is comparable to our fastest de- 
clining sub-Chandrasekhar-mass model and also too fast 
compared to normal SNe la. 

In addition, we note that Am\§{E) is also quite sen- 
sitive to the details of the radiative transfer treatment. 
Using an ato mic data set with 8 x 10 6 lines (the big_gf- 
4 data set of IKromer fc Simll2009D . for example, yields 
Am\5{B) ~ 1.75 for W7. In contrast, a simulation using 
the atomic data set adopted in this study but applying a 
pure LTE treatment for the excitation/ionization state of 
the plasma yields Am\${B) ~ 1.95 for W7. This shows 
that systematic uncertainties in the radiative transfer 
treatment can affect Amis(i3) by several tenths of a 
magnitude (see also the compari son of different radia - 
tive transfer codes in Figure 7 of IKromer fc"Siir] [2009ft . 
Thus, we argue that there is no evidence that our sub- 
Chandrasekhar-mass models fade too fast compared to 
other explosion models. 

4.2. Colour evolution 

The most striking difference between our light curves 
and those of the comparison objects in Figure [5] con- 
cerns their colour. To highlight this we show the angle- 
averaged time-evolution of the B — V, V — R and V — I 
colours for all our models in Figure |3] and compare them 
again to our fiducial SNe 2005cf, 2004eo, 2001el and 
1991bg. 

In B — V all our models show positive colour indices 
for the whole simulation period and a red peak at ~ 40 
days after explosion. Contrary to observed SNe la ([Lira! 
1995), however, we find no convergence of the different 
models at epochs after the red peak. Instead, our mod- 
els at all times form a sequence of increasingly redder 
B — V colour towards the fainter explosions. With the 
exception of Model 6, all our models are generally too 
red compared to spectroscopically normal SNe la. At 
maximum light, for example, we find {B — U)-values of 
1.67 and 0.28 for Model 1 and Model 5, respectively (val- 
ues for the other models are given in Table [T]) . In con- 
trast, spectroscopically normal SNe la are characterized 
by (B - V) max ~ 0.0 mag. 

Subluminous 1991bg-like objects show a redder B — 
V colour before the red peak, r eaching B — V ~ 
0.4 . . . 0.7 mag at B-band maximum ([Taubenberger et al.l 
2008). Although Model 4 can reproduce this, it is not a 
good fit to 1991bg-like objects, since it is considerably 
too bright (compare Figure [2]). Our subluminous Models 
1 and 2, on the other hand, are significantly redder than 
observed {B — V colours of 1.67 and 1.51 at maximum , 
resp ectively) . In contr ast, IHoflich fc Khokhlovj (|1996D 
and IHoflich et al.l (|1996[ ) found too blue colours at max- 
imum light for their subluminous sub-Chandrasekhar- 
mass models compared to 1991bg-like objects. 

The origin of our red colours traces back to the ex- 
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Figure 4. Angle-averaged (thick black line) spectra at three days before £?-band maximum for all six of our models as indic ated by 
the labels. For comparison the blue line shows the de-redshifted and de-reddened spectrum of SN 2004eo (Pastorcllo ct al. 2007a) at the 
corresponding epoch. This was scaled such that its maximum matches the maximum of the model spectrum. The colour coding indicates 
the element(s) responsible for both bound-bound emission and absorption of quanta in the Monte Carlo simulation. The region below the 
synthetic spectrum is colour coded to indicate the fraction of escaping quanta in each wavelength bin which last interacted with a particular 
element (the associated atomic numbers are illustrated in the colour bar). Similarly, the coloured regions above the spectra indicate which 
elements were last responsible for removing quanta from the wavelength bin (either by absorption or scattering/fluorescence). 



tended layer of titanium and chromium which is present 
in the helium shell ejecta of our models. This will be 
discussed in more detail in Section l4~3l 

While the trend of increasingly redder colours for 
fainter models persists in the V — R index, our models are 
not systematically too red here. Model 2 and 3 populate 
about the right range for 1991bg-like and spectroscopi- 
cally normal SNe la, respectively. However, the details 
of their V — R evolution, especially the initial decline 
present in all the models, do not match the observations. 

A similar behaviour is found for the V ~I colour where 
again Model 2 and 3 lie closest to the observed colours 
of 1991bg-like and spectroscopically normal SNe la, re- 
spectively. However, again the agreement is imperfect. 
For Model 2 this is most obvious at the latest epochs, 
where 1991bg-like objects start to become bluer while 
the model colour maintains a redward evolution. More- 
over, the model shows a secondary blue minimum, due 
to the post-maximum plateau of the T^-band light curve 
which is not observed in 1991bg-like objects. In Model 3 
the initial decline and the rise to the red peak are signif- 
icantly different from the observed behaviour of spectro- 
scopically normal SNe la. 

4.3. Spectra 

The properties of our model light curves can be un- 
derstood from consideration of the synthetic spectra. 
Figure 0] shows the angle-averaged spectrum of all our 
models at three days before -B-band maximum. For 
comparison, we also show the spectrum of SN 2004eo 



(jPastorello et al.l l2007af) at the corresponding epoch. 
The colour coding below the synthetic spectrum indi- 
cates the fraction of escaping quanta in each wavelength 
bin which were last emitted in bound-bound transitions 
of a particular element. Similarly, the coloured regions 
above the spectra indicate which elements were last re- 
sponsible for removing quanta from a wavelength bin by 
bound-bound absorption. This coding allows us to both 
identify individual spectral features and track the effect 
of fluorescence on the spectrum formation directly. The 
contributions of bound-free and free- free emissions to the 
total flux are so small that they are not discernable in 
the figure. 

From this it is immediately obvious that the colours 
of our models are due to fluorescence in titanium and 
chromium lines. Having a wealth of strong lines in the 
UV and blue, even small amounts of these elements ef- 
fectively block radiation in the UV and blue and redis- 
tribute it to redder wavelengths where the optical depths 
are smaller and the radiation can escape. Compared to 
other explosion models, this effect is particularly strong 
in the lFink et alJ J2010) models, since they produce rel- 
atively large amounts of titanium and chromium in the 
outer layers during the initial helium detonation (cf. Ta- 
bleland Figure [TJ typical titanium and chromium yields 
for other explosion models are on the order of the yields 
from the core detonation). This also explains the trend 
for redder colours in the fainter models: according to Ta- 
ble [TJ the production of titanium and chromium in the 
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shell increases continuously from Model 6 to Model 1. 

Since this titanium and chromium layer is located at 
higher velocities than most of the intermediate mass 
elements (cf. Figure Q}, redistribution by titanium 
and chromium also dilutes the absorption features of 
intermediate-mass elements like silicon and sulphur by 
reprocessing flux into the relevant wavelength regions. 
This can be clearly seen in Figure |4j Although Models 
1 and 2 produce the largest amounts of silicon of all our 
models (0.27 and 0.25 solar masses in the core, respec- 
tively), they show only a (very) weak Sin A6355 feature 
since this wavelength region is strongly polluted by flux 
redistribution from titanium and chromium. The only 
feature of intermediate-mass elements which is clearly 
visible in these models is the Can near-infrared (NIR) 
triplet at AA8498, 8542, 8662. This feature remains un- 
buried since calcium exists co-spatially with titanium and 
chromium in the outer shell (Figure [TJ. 

Of all our models, Model 3 [M core ( 28 Si) = 0.21 M Q ] 
produces the strongest Sin A6355 feature. For more 
massive models the silicon yields from the core detona- 
tion drop dramatically since, owing to the higher den- 
sities, a larger fraction of the core is burn ed to iron- 
group elements (see Table [1] and Figure 3 of iFink et al.l 
120101 ). As expected from st udies of pure deton ations of 
Chandrasekhar-mass WDs ([Arnett et al.lll97lh . the ex- 
treme case of Model 6 (almost the Chandrasekhar mass) 
produces no significant amounts of intermediate-mass el- 
ements. Thus the spectrum is totally dominated by iron- 
group elements. Since it shows no indication of any Si n 
A6355 feature, this model would not be classified as a 
SN la. 

Since none of our models give a good match to the 
colours and line strengths of observed SNe la, we will not 
discuss line velocities in detail. We note, however, that 
along the sequence from Model 1 to 5 there is a trend for 
higher line velocities of intermediate mass elements. This 
is obvious in the Sin A6355 feature of our models which 
is apparent in Figure |U Along the model sequence (1 
to 5), this line moves to shorter wavelengths compared 
to the observed Sin absorption feature of SN 2004eo. 
This arises because the inner boundary of the region rich 
in intermediate-mass elements moves to higher velocities 
with increasing mass due to the more c omplete burning i n 
the inner regions (compare Figure 3 of lFink et alJl2010h . 

Finally, we note that none of our model spectra show 
any indication of helium lines, despite our models having 
up to ~ 0.08 Mq of helium in their outer layers (cf. Ta- 
bled]). Since helium has rather highly excited levels, this 
might simply be a consequence of our approximate treat- 
ment of the plasma state which neglects non-thermal ex- 
citation and ionization. We note, however, that despite 
using the phoenix code which does include a treatment 
of non-thermal processes. iNugent et a l. (1997) also found 
no strong evidence of helium lines for models with even 
more (0.2 Mq) helium in the outer shell. 

4.4. Line-of-sight dependence 

As discussed in Section [2j our models show strong 
ejecta asymmetries due to their ignition in a single-point 
at the north pole of the WD (cf. Figure [T]). These asym- 
metries are expected to have some influence on the ob- 
servables along different lines-of-sight. Since they are 
characteristically the same for all our models, we first 



1 1 






cos 9 


-o.9 : 

0.0 1 
0.9 


lr" 


if 
1 vJ 



















2000 4000 6000 8000 

Wavelength in A 



Figure 5. Line-of-sight dependent maximum-light spectra for 
Model 3. To indicate the maximal effect spectra seen pole-on are 
plotted and compared to a spectrum seen equator on. The corre- 
sponding lines are identified by the colour coding. For comparison 
the angle-averaged spectrum is shown as the dashed black line. 

use Model 3 as an example to give a detailed discus- 
sion of line-of-sight dependent spectra and light curves 
of this model fSections f4.4.1l and l 4.4.2j) . In Section l4~4~3l 
we then present the variations in peak magnitudes and 
colours due to line-of-sight effects for all our models. 

4.4.1. Spectra 

To obtain line-of-sight dependent observables we bin 
the escaping photons into a grid of ten equal solid-angle 
bins in /i = cos with 6 being the angle between the line- 
of-sight and the z-axis of the model. In Figure [5] we show 
the maximum-light spectra of Model 3 for three differ- 
ent directions, corresponding to lines-of-sight close to the 
southern polar axis (/i = —0.9), the northern polar axis 
(fi = —0.9) and equator-on (average of the fi = —0.1 and 
0.1 bins). While the equator-on spectrum looks similar 
to the angle-averaged spectrum, the spectra seen from 
the polar directions are significantly different from each 
other as well as from the angle-averaged spectrum. For 
/i = 0.9, an observer is looking through the extended 
layer of iron-group elements in the outer shell of the 
northern- hemisphere (cf. Figure [T]) . This blocks almost 
all the flux in the UV and blue wavelength range and 
the redistribution of this flux leads to a very red spec- 
trum. In contrast, this layer of iron-group elements in 
the outer shell is far less extended on the southern hemi- 
sphere, where the shell burning is less efficient. Thus, 
an observer with [i — —0.9 sees a bluer spectrum than 
equator-on. Equator-on the extension of the layer of iron- 
group elements is somewhere between these two extremes 
and the spectrum resembles the angle-averaged case. 

A secondary effect results from the off-centre ignition 
of the core. Since the core has been compressed less 
strongly on the northern hemisphere, the core detonation 
yields more intermediate-mass elements on the northern 
than on the southern hemi sphere (cf. Figu re [1] and the 
Discussion in Section 4.3 of lFink et alll2010h . This leads 
to stronger features of intermediate-mass elements in the 
spectrum seen along the northern polar-axis and becomes 
particularly visible in the strength of the Si n 6355 A fea- 
ture and the Ca n NIR-triplet which become weaker from 
^ = 0.9^-0.9. 
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Figure 6. Selected line-of-sight dependent light curves of Model 3 as indicated by the colour coding. For comparison angle- averaged light 
curves (black dashed) and photometrical data of our fiducial SNe 2005cf, 2004eo, 2001el and 1991bg (different symbols) are shown. 



4.4.2. Lightcurves 

Figure [5] shows band-limited synthetic light curves 
from our radiative transfer simulations for Model 3 as 
seen equator-on and from the two polar directions. As 
already noted for the maximum-light spectra, our model 
shows a trend for increasingly red colours from /i = 
—0.9 — > 0.9. Similarly, we find a clear dependence of 
the light curve rise and decline times on the line-of- 
sight. While the rise times increase along the sequence 
[i = —0.9 — > 0.9 (for B band, for example, we find 
rise times between 16.8 and 19.7 days for fi = —0.9 
and fi = 0.9, respectively), the light curve declines 
more slowly [Ami5(B) = 1.91 and /S.m\^{B) = 1.29 for 
/i = —0.9 and /i = 0.9, respectively]. Moreover, we find 
a clear trend for a weaker dependence of the light curves 
on the line-of-sight at lower photon energies. Thus, while 
U and B band show a variation of ~ 3 and ~ 2 magni- 
tudes at maximum light, respectively, the variation in V 
band is already less than half a magnitude and in the 
NIR bands there is virtually no viewing-angle effect. 

All these effects are due to the asymmetry of the layer 
of iron-group elements in the outer shell. Since this layer 
is more extended on the northern hemisphere it causes 
additional line blocking and thus enhanced fluorescence 



and photon trapping for inclinations close to \x = 1. This 
explains the redder colours as well as the increasing rise 
and decreasing decline times for fi = —0.9 — > 0.9. For 
lower photon energies, the asymmetry of the outer shell 
is less important, since the optical depths are smaller 
and photons typically escape from deeper layers of the 
ejecta. In the NIR bands the entire ejecta contribute to 
the emitted photons and the viewing-angle dependence 
of the light curves is very small. 

This is illustrated in Figure [7] which shows where pho- 
tons of selected bands were last emitted before they es- 
caped from the supernova ejecta at different times. While 
[/-band photons leak out predominantly from the re- 
gions on the southern hemisphere, where the layer of 
iron-group elements is least extended, /-band photons 
show no strong preference for a particular region. In 
fact, even before maximum light the whole ejecta con- 
tribute to /-band emission. V^-band photons, in contrast, 
show some preference for leaking out from the southern 
hemisphere before and immediately after maximum light. 
From about 10 days after maximum light this preference 
disappears. This is directly reflected in the y-band light 
curve in Figure |6] which becomes viewing-angle indepen- 
dent at about that time. 
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Figure 7. Region of last emission for selected bands (U, V, I from top to bottom) and different times (from left to right). The epochs 
indicated in the bottom panels are given in days with respect to B-band maximum. Dark regions contribute most to the flux escaping in 
the band. The model is symmetric under rotation about the z-axis. 



Finally, we note that the strong line-of-sight depen- 
dence of our sy nthetic observabl es poses an additional 
problem for the iFink et al.l (|2010| ) models. Even if some 
particular line-of-sight might compare more favourably 
to observed SNe la than others (as does for example the 
southern line-of-sight for the light curves in Figure [B]), all 
lines-of-sight should occur in nature - including the most 
peculiar ones. Moreover, the large dispersion in bright- 
ness at -B-band maximum which we find in our model is 
in conflict with observations. 

4.4.3. Other models 

Since the other models in the IFink et ahl (|2010T l se- 
quence have the same characteristic asymmetries, their 
light curves and spectra show a similar viewing-angle de- 
pendence. However, the strength of this viewing-angle 
dependence varies between the models, due to their dif- 



ferent helium shell masses. To demonstrate this, Figure[8] 
shows selected light curve properties of all models for the 
ten different viewing directions within our uniform grid 
in fi. 

As expected, the viewing-angle dependence decreases 
for smaller helium shell mass. Thus, we find a scatter 
of more than 2 magnitudes for the brightness at B-band 
maximum between the different lines-of-sight in Model 
1 and 2, while Model 6 shows only a scatter of ~ 0.3 
magnitudes. Similar trends can be observed for the B — V 
and V — R colours at £>-band maximum. While Model 
6 shows almost no dependence on the line-of-sight, the 
less massive models show a very clear trend for redder 
colours towards /i = 0.9 (as discussed for Model 3 in 
detail above). 

For the rise times the situation is similar: while Models 
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Figure 8. Light curve properties of the Fink ct al. (2010) models 
for 10 different viewing directions (different models are indicated 
by the colour coding). From top to bottom the panels show the 
rise-time to S-band maximum [t ma x(-B)], the decline in B-band 
between maximum light and 15 days thereafter [Ami5(B)], the 
B-band peak magnitude (Mb, max), the B — V colour at B-band 
maximum [(B-V)B max ] an d the V — R colour at S-band maximum 
[(V — i?)Bmax]> respectively. The lines connecting the different data 
points are just to guide the eye. 



1, 2 and 3 (relatively massive shells) show a clear trend 
for increasing i3-band rise times from fi = —0.9 to fi = 
0.9 due to the enhanced photon trapping by optically 
thick lines of iron-group elements, we do not observe a 
clear trend for Models 4, 5, and 6 (less massive shells). 
Recall, that the average rise times increase from Model 1 
to 3, while Models 4 to 6 form a sequence of decreasing 
rise times. This was already discussed in Section 03] 

The enhanced trapping of photons by the iron-group 
layer of the northern hemisphere also affects the post- 
maximum decline rate in B band. Thus, we find decreas- 
ing values of Amis (I?) from \i = —0.9 to /i = 0.9. This 
trend persists for all our models, although it is weaker 
for the models with the least massive shells. 




20 40 60 

Time since explosion (days) 




0.1 1.0 

Photon energy (MeV) 

Figure 9. 7-ray light curves (top panel) and sp ectra (lower pan- 
els) for Models 1, 3 and 6 of Fin k et al. (2010) as indicated by 
the colour coding. The spectra shown are for two different epochs 
which correspond roughly to maximum light in B band (20 days, 
middle panel) and in 7-rays (60 days, bottom panel). Line features 
discussed in the text are identified by labels (A - F) in the middle 
panel (for their identification, see Table I2J. In addition, the 56 Ni, 
56 Co and 48 V line systems are indicated. The line systems of 48 Cr, 
52 Fe and 52 Mn are omitted, since their life times are sufficiently 
short that these nuclei have already decayed at the shown epochs. 
Note that at 20 days the line features in the spectra are generally 
offset towards higher energies compared to the line identifications, 
due to blue-shifted emission. At 60 days the emission comes from 
regions of lower velocity and the offset disappears. 



4.5. 7-ray emission 

Due to their peculiar composition including a mixture 
of the radioactive isotopes 56 Ni, 52 Fe and 48 Cr close to 
the surface, 7-observations might provide an additional 
discriminant between our sub-Chandrasekhar-mass mod- 
els and more standard explosion models which do not 
show radioactive isotopes close to the surface. To inves- 
tigate this, Figure |9] shows 7-ray light curves and spectra 
for Models 1, 3 and 6 of our sequence. Line identifica- 
tions are given for some important features in Table [5] 

Broadly speaking, our 7-ray spectra a re not dramat- 
ically different from those obtained by iSim Sz Mazzalil 
(2008) for paramcterised Chandrasekhar-mass models: 
they are dominated by strong emission lines, mainly due 
to Co, and a continuum which results from Compton 
scattering of line photons. 

Models 1 and 3 have characteristically similar 7-ray 
light curves. After an initially fast rise (lasting for about 
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Table 2 

Identification of 7-lines in Figure [9] 



Identifier Photon energy (MeV) Source 



A 0.158 56 Ni 

B 0.270 56 Ni 

C 0.511 Annihilation of positrons from 56 Co and 48 V 

D 0.750 56 Ni 

E 0.812 56 Ni 

F 1.562 56 Ni 



5 and 10 days after the explosion for Model 3 and 1, re- 
spectively) their light curves have a small plateau before 
passing to a second rise to maximum light at about 60 
days. In contrast, Model 6 shows no early rise/plateau 
but reaches maximum at about 60 days in one continuous 
rise. These differences result from the significantly dif- 
ferent masses of the helium shells of the models. Model 1 
has a rather massive shell of 0.126 M© with ~ 0.02 M© 
of radioactive isotopes. The fast initial rise of the light 
curve of Model 1 is caused by 7-photons which originate 
from this outer shell. Due to the small optical depth in 
the outer shell, those 7-photons start to stream freely at 
about 10 days and the emerging flux from the outer shell 
decreases. At the same time, however, 7-photons from 
the C/O core start to escape and keep the light curve ris- 
ing until the core also becomes transparent to 7-rays at 
about 60 days after which the 7-photons start to stream 
freely. 

For Model 3, we observe the same effect. However, 
due to its lower helium shell mass (M s h = 0.055 M©) 
7-photons from the shell escape even earlier (at about 
5 days). Since the detonation of the C/O core of this 
model produces much more 56 Ni, the 56 Ni-rich region is 
far more extended than in Model 1. This means that 56 Ni 
nuclei from the C/O core are present in regions where the 
Compton optical depth is relatively low, allowing their 
7-photons to more easily escape. Thus, photons from the 
core dominate the rise earlier than in Model 1. In the 
extreme limit of Model 6, the shell (M sh = 0.0035 M ) 
is completely negligible and the massive 56 Ni region in 
the C/O core extends so close to the surface that the 
rise of the light curve is totally dominated by 7-photons 
escaping from the C/O core. 

A similar effect shows up in the early-time (20 days) 
7-ray spectra of our models. While Model 6 shows a 
clear indication of the 0.158 and 0.270 MeV lines of 
56 Ni, those are invisible for Model 1. Since the Comp- 
ton cross-section increases with decreasing photon en- 
ergy, soft-energy lines are most easily buried by photons 
being Compton down-scattered from higher energies and 
can onl y be observed if 56 Ni is p resent at low optical 
depths |G omez- Gomar et aLlll998D . Thus the presence 
of the 0.158 and 0.270 MeV lines of 56 Ni in Model 6 is 
a direct consequence of the large extension of the 56 Ni 
bubble in this model, compared to Model 1. The harder 
56 Ni lines at 0.750, 0.812 and 1.562 MeV, in contrast, are 
also visible in Model 1. However, they are weaker due to 
the smaller mass of 56 Ni synthesized in Model 1. Note 
that despite containing 48 Cr or 48 V close to the surface, 
our models show no clear features of these radioactive 
isotopes in their 7-ray spectra. 



An interesting effect concerns the 511 keV annihilation 
line: for Model 3 and 6 the strength of this line increases 
from 20 to 60 days significantly, but it does not for Model 
1. In Model 1, the 511 keV line is dominated by positrons 
from 48 V at 20 days. Being located in the outer layers, 
the annihilation photons escape easily making the line 
strong. At 60 days most of the 48 V has already decayed. 
Then, the 511 keV line results from the annihilation of 
56 Co positrons from the C/O core. For Model 3 and 6, 
in contrast, the line is always dominated by annihilation 
photons originating from 56 Co positrons in the C/O core. 
Due to the longer life time of 56 Co and the longer time it 
takes for photons to escape from the core, the strength 
of the 511 keV line increases from 20 to 60 days. 

5. DISCUSSION 

In the last Section we presented synthetic observables 
for the sub-Chandrasekhar-mass doub le detonation mod - 
els with minimum helium shell mass of lFink et al.l (2010). 
Compared to observed SNe la, these models show some 
promising features: 

1. They predict a wide range of brightnesses that cov- 
ers the whole range of observed SNe la. 

2. Their light curve rise and decline times are for 
most bands in reasonable agreement with those ob- 
served. 

But despite these positive features the iFink et aLl 
(2010) models cannot account for all the properties of 
observed SNe la since they have peculiar light curves 
and spectra: 

1. The colours of the models are too red compared to 
observed SNe la. This is particularly obvious in the 
evolution of the B — V colour, where all the models 
are redder than spectroscopically normal SNe la for 
all epochs. 

2. The model spectra cannot reproduce the strong 
features of intermediate-mass elements typical of 
SNe la at maximum light. 

In detail, there are further problems concerning the 
exact light curve shapes and decline rates as well as a 
strong viewing-angle dependence which is caused by the 
point-like ignition of our models. 

We have argued that all these problems are mainly 
due t o the burn i ng pr oducts of the helium shell. More- 
over, iSim et akl (|20100 . have recently shown that deto- 
nations of centrally ignited spherically symmetric naked 
sub-Chandrasekhar-mass C/O WDs are capable of repro- 
ducing the observed diversity of light curves and spectra 
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Figure 10. Selected line-of-sight dependent light curves of Model 3c as indicated by the colour coding. For comparison angle-averaged 
light curves (black dashed) and photometrical data of our fiducial SNe 2005cf, 2004eo, 2001el and 1991bg (different symbols) are shown. 



of SNe la - at least to a similar level of agreement as 
models of the more standard Chandrasekhar-mass de- 
layed detonations. This naturally leads us to speculate 
on whether the helium shell in double detonation mod- 
els could be altered in some way to reduce its negative 
impact on the synthetic observables. 

5.1. Influence of the helium shell 

We first explicitly investigate the extent to which the 
shortcomings and viewing-angle dependence of our mod- 
els can indeed be attributed to the helium shell. For that 
purpose we constructed (for Model 3) a toy model which 
contains only the burning products of the detonation in 
the initial C/O core but not those o f the i nitial helium 
shell. For the models of iFink et al.l (|2010() this can be 
done in a straightforward manner: since the models use 
two different sets of tracer particles to simulate the nu- 
cleosynthetic yields of core and shell b urning respectively 
(see Section 3.3 of IFink et all [2010). we obtain such a 
model by restricting our SPH-like reconstruction algo- 
rithm (Section [3]) to the core tracers. Properties of this 
"core-only" model (hereafter 3c) are listed in Table [1] 

Figure [10] shows band-limited synthetic light curves ob- 
tained from our radiative transfer simulations for this 



model as seen equator-on and from the two polar di- 
rections. In contrast to the light curves of Model 3 in 
Figure [6[ which are strongly dependent on the viewing 
angle, the light curves of Model 3c show only a moder- 
ate line-of-sight dependence. Thus, at maximum light 
we find now only a variation of ~ 0.5 and ~ 0.2 mag- 
nitudes for U and B band, respectively. For Model 3 
these values were significantly larger (~ 3 and ~ 2 mag- 
nitudes, respectively). Redder bands show no significant 
line-of-sight dependence in Model 3c. 

Moreover, the light curves of Model 3c give an excellent 
representation of SN 2004eo in the B, V and R bands. U 
and I are not in perfect agreement, but still reasonable 
compared to the agreement between other fi rst principles 
explo sion models and observed SNe la (e.g. iKasen et al.l 
2009). In particular, the colours of this toy model are 
now in good agreement with observed SNe la and not 
too red as it is the case for Model 3. In the NIR bands, 
in contrast, the agreement is no better. However, the 
NIR light curves are much more difficult to model ac- 
curately since they require simulations with an exten- 
sive atomic data set to properly simulate flux redistribu- 
tion by fluorescence which st rongl y affects these bands 
(jKasenl [20061 : iKromer fc Siml 12009( 1 . Here, however, we 
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Figure 11. Angle-averaged (thick black line) spectra at three 
days before £?-band maximum for Model 3c. For comparison the 
blue line sh ows the de-redshifted and de-reddened spectrum of 
SN 2004eo (Pastorcllo ct al. 2007a) at the corresponding epoch. 
Note, that the flux is here in physical units and not scaled like in 
Figure [4] For a description of the colour coding see Figure [4] 

have restricted ourselves for computatio nal reasons to a 
simp lified atomic data set (cd23_gf-5 of iKromer fc Siml 
120091 ) with only ~ 400,000 lines. This has been shown 
to give reliable results in the optical bands which, for 
our purposes, are the most important since they are the 
most differ e nt be tween our toy Model 3c and Model 3 of 
iFink et al.l (|2010D . In the NIR bands, in contrast, Model 
3c and Model 3 give rather similar results (see Figure [T2)) . 

The good agreement between our toy Model 3c and 
observational data does not only hold for band-limited 
light curves but also for individual spectral features as 
can be seen from Figure [TT] which shows a spectrum of 
Model 3c at 3 days before B-band maximum. Compared 
to SN 2004eo, our toy model succeeds in reproducing the 
characteristic spectral features of intermediate-mass ele- 
ments in SNe la. This is highlighted by our colour coding. 
Moreover, it shows an overall flux distribution which is 
in almost perfect agreement with the observed spectrum 
of SN 2004eo and we see no strong flux redistribution by 
titanium (compare with Model 3 in Figure [4]) . 

This confirms our conclusion from Section l4~4l that the 
peculiarities of our model spectra with respect to the 
observations and their strong viewing-angle dependence 
are mainly due to the shell material and its compositional 
asymmetries. It also shows that the off-centre ignition of 
the secondary detonation in the C/O core causes only a 
minor viewing-angle dependence which is on the order of 
the observed variation between SNe la. 

5.2. Prospects 

In light of the discussion above, we are motivated to 
speculate on how the influence of the helium shell might 
be reduced. In the sub-Chandrasekhar-mass double det- 
onation scenario, the helium shell cannot be removed en- 
tirely since it is required to trigger th e detonation. Also 
the helium shell mass adopted in the IFink et ah! ([2010D 
models is alr eady the minimum th at might be expected 
to detonate (jBiidsten et al.l 120071) . In Section H how- 
ever, we have argued that the differences between our 



model spectra and observations are not a consequence of 
the helium itself but of its particular burning products, 
namely titanium and chromium produced in the outer 
layers. The yields of these elements are affected by de- 
tails of the nucleosynthesis in the shell. 

The degree of burning in the shell material (and thus 
its final composition) can be affected by the initial 
abundance of heavy nuclei (e.g. 12 C) which in turn de- 
pend strongly on triple-a reactions during previous hy- 
drostatic burning and dredge-up phases from the core 
(|Shen fc Bildstenl [2009). Since the time-scale for a- 
captures behind the detonation shock front is signifi- 
cantly shorter than that of triple-a reactions, such seed- 
nuclei can limit the a-chain before reaching nuclear sta- 
tistical equilibrium. If, for example, in a shell consisting 
of a mixture of He and 12 C the number ratio of free 
a-particles to 12 C-nuclei on average is less than 6 (cor- 
responding to a mass ratio of 2), the a-chain will end 
at 36 Ar. Thus, it is possible that more intermediate- 
mass elements and less titanium and chromium may be 
produced. Therefore it is interesting to consider how 
the burn ing of the helium might be different from that 
found by IFink et al.l (|2010D for different initial composi- 
tions of the shell. A full study of this goes beyond the 
scope of this work and will be published in a follow-up 
study. Here, we illustrate the possibility of obtaining 
better ag reement with data for just one example. 

In the IFink et ail (|2010l ) models it was initially as- 
sumed that the shell consisted of pure helium. To demon- 
strate the sensitivity to the initial composition of the 
shell, we set up another toy model. For this "modified" 
model (hereafter 3m), we homogeneously polluted the 
shell of Model 3 with 34% (by mass) of 12 C and repeated 
the hydrodynamics and nucleosynthesis calculation (in 
the same way as described by IFink et aLI [20100 We 
found that a core detonation was still triggered but the 
different shell burning led to a substantial reduction of 
the mass of 44 Ti, 48 Cr and 52 Fe in the shell (nucleosyn- 
thetic yields for core and shell of the mode l are given 
in Tab le [J) . Since the detonation tables of Fin k et al.l 
( 2010) are only valid for pure helium, Model 3m is not 
fully self-consistent. Nevertheless, it is a useful toy model 
to explore the basic effect of a modified shell composition. 

Figure [12] compares the angle-averaged band-limited 
light curv es of this modifie d model (3m), to those of 
Model 3 of lFink et al.l (|2010D and our core-only toy model 
(3c). As can be seen, the modified Model 3m produces 
light curves very similar to those of Model 3c despite 
having about the same shell mass as Model 3. The most 
obvious difference between our modified and core-only 
models occurs in the U band: the titanium in the outer 
layer of Model 3m causes some line blocking, leading gen- 
erally to a dimmer £/-band magnitude than for Model 3c 
which has no outer layer. Compared to Model 3, with 
its large titanium mass in the shell, however, this effect 
is much weaker. The B band, which is strongly affected 
in Model 3, shows no significant titanium absorption for 
Model 3m. Another slight difference between Model 3m 

1 We note, that the base temperature of the shell in Model 3m 
was decreased to 4 X 10 s K (compared to 6.7 X 10 s K in Model 3) 
to suppress further triple-a burning in the shell. As a consequence 
the density and also the shell mass changes slightly compared to 
the original Model 3 (cf. Table [TJ. 
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Figure 12. Angle-averaged UVOIR bolometric and U ,B,V ,R,I ,J ,H ,K band-limited light curves for Model 3, 3c and 3m of our model 
sequence as indicated by the colour coding (compare Table [T] for details on the models). For comparison angle-averaged light curves (black 
dashed) and photometrical data of our fiducial SNe 2005cf, 2004eo, 2001el and 1991bg (different symbols) are shown. 

and 3c occurs after the first peak in the I band. Com- 
paring the light curves of Model 3m to SN 2004eo, we 
find qualitatively similar agreement as for Model 3c. 

This generally also holds for the angle-averaged spec- 
trum at 3 days before £?-band maximum, shown in Fig- 
ure [13l Compared to Model 3c, where the agreement 
was almost perfect, there are some minor shortcomings. 
But the m odel is dramatica lly improved compared to 
Model 3 of iFink et all (|2010D . The small differences be- 
tween Model 3c and 3m are again mostly due to the 
titanium in the outer layers which leads to pronounced 
absorption troughs bluewards of the Ca II H and K lines 
and redwards of 4, 000 A. This suggests that Model 3m 
still over-produced titanium in the shell. Interestingly, 
the enhanced calcium abundance in the outer layers (cf. 
Table [IJ leads to a stronger Can NIR triplet, bring- 
ing the model in better agreement with the spectrum of 
SN 2004eo at the corresponding epoch. Therefore some 
calcium in the outer shell is an improvement over Model 
3c. This suggests that a slight further reduction in the 
degree of burning so that titanium is further suppressed 
in favour of calcium (just one step down the a-chain) 
could lead to very good agreement. 
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Figure 13. Angle- averaged (thick black line) spectra at three 
days before B-band maximum for Model 3m. For comparison 
the blue line shows the de-redshifted and de-reddened spectrum 
of SN 2004co (Pastorcllo ct al. 2007a) at the corresponding epoch. 
Note, that the flux is here in physical units and not scaled like in 
Figure [4] For a description of the colour coding see Figure HI 
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In summary, polluting the initial helium shell of Model 
3m with 12 C significantly improved the agreement be- 
tween our synthetic spectra and light curves and those 
observed for SNe la making this model a promising pro- 
genitor candidate for SNe la. We stress again that this 
improvement results only from the change in the com- 
position of the burning products of the helium shell 
which contains much less titanium and chromium for 
this model (the total shell mass stays about the same). 
Given that the initial composition of the helium shell 
depends on several processes including details of the 
accretion physics and hydrostatic burning phases that 
might precede the detonation or possible dredge-up of 
material from the C/O core, this leaves some scope to 
find sub-Chandrasekhar-mass double detonation models 
which give reasonable agreement with observed SNe la. 
This, however, must be investigated by future follow-up 
studies that more fully explore the influence of the initial 
composition of the helium shell on the burning products 
and link the initial composition of the helium shell di- 
rectly to the evolution of progenitor models. Moreover, 
we note that different ignition geometries might also lead 
to better agreement with observational data. In partic- 
ular, more symmetric ignition geometries, e.g. ignition 
in an equatorial ring or simultaneous ign ition in multi- 
ple points (as studied bv lFink et al.ll2007l for the case of 
more massive helium shells), are likely to alleviate the 
strong viewing -angle dependence found for the point-like 
ignition of the lFink et al.1 {2010) models. 

Our results also highlight the strong sensitivity of the 
radiative transfer to particular elements/ions (in our case 
titanium and chromium which represent only a tiny frac- 
tion of the ejecta mass yet dominate our conclusions). 
This emphasizes the need for a better description of nu- 
clear reaction rates and continued study of the radiative 
transfer processes (and atomic data) in order to quantify 
more fully the systematic uncertainties which arise due to 
the complexity of spectrum formation in supernovae. In 
particular, we note that almost all the flux redistribution 
done by titanium and chromium in our models is due to 
their singly ionized states. Since the current ionization 
treatment of ARTIS neglects non-thermal processes (see 
IKromer fc Shr][20Q9l for more details), we cannot exclude 
that the actual ionization state in the helium shell ejecta 
would be higher due to non-thermal ionization from the 
radioactive isotopes produced during the helium burning. 
This could also significantly improve the agreement be- 
tween our models and observational data, as a numerical 
experiment with an artificially enforced higher ionization 
state for titanium and chromium has shown. 

6. CONCLUSION 

In this paper we presented synthetic observables for 
th e sub-Chandrasek har-mass double detonation models 
of iFink et all (|2010() . We found that these models pre- 
dict light curves which rise and fade on time-scales typ- 
ical of SNe la. Moreover, they produce a large range of 
brightnesses which covers the whole range of observed 
SNe la. However, they do not account for all the proper- 
ties of observed SNe la since they have peculiar spectra 
and light curves. In particular, their B — V colours are 
generally too red compared to observed SNe la. This is 
in contrast to the results of e arlier work on models with 
more massive helium shells (IHoflich fc Khokhlovl 11996c 



IHoflich et aHl!996t iNugent et alj|1997l ). In addition, our 
model light curves and spectra show an unreasonably 
strong viewi ng-angle dependen ce due to the point-like ig- 
nition of the IFink et aLl (|2010| ) models and the resulting 
ejecta asymmetries. 

Detonation of a pure helium shell leads to a layer con- 
taining iron-group elements like titanium and chromium 
around the core ejecta. These elements have a vast num- 
ber of optically thick lines in the UV and blue part of 
the spectrum making them very effective in blocking the 
flux in these wavelength regions and redistributing it to 
the red. We used a toy model to show that this layer of 
titanium and chromium causes the peculiar red colours 
of our light curves and also the peculiar spectral fea- 
tures. Moreover, we found that this toy model repro- 
duces the observed properties of SNe la remarkably well. 
The toy model also showed that the strong viewing-angle 
dependence of our models results from the compositional 
asymmetry in the helium shell ejecta and not from the 
off-centre ignition of the C/O core. We stress, that the 
additional energy release in the shell, due to the produc- 
tion of radioactive nuclides during the helium burning, is 
relatively inconsequential for our models - even at 7-ray 
energies the signatures of the surface 48 Cr and 52 Fe are 
not apparent. Instead, in the optical/UV the shell has 
a strong signature but this is primarily due to the addi- 
tional opacity in the outer layers which affects the trans- 
port of energy from the core to the surface. We conclude 
that, if the double detonation sub-Chandrasekhar-mass 
model valid for normal SNe la, the properties of the post- 
burning heli um shell material need to be different from 
those in thelFink et all (120101) models. 

bmce IFink et al.l ifeoiO) considered the limit of the 
least massive helium shells which might ignite a deto- 
nation in the helium shell, their models represent the 
most optimistic case for reducing the influence of the 
shell material by simply reducing the shell mass. How- 
ever, we argue that the mass of the helium shell ejecta 
is not the main problem but rather the peculiar com- 
position including comparably large masses of titanium 
and chromium. We illustrated this using a second toy 
model, where the initial composition of the helium shell 
was polluted with 34% (by mass) of 12 C. By providing ad- 
ditional seed-nuclei for a-captures, this leads to burning 
products with lower atomic number (i.e. intermediate- 
mass elements rather than mainly iron-group elements). 
Spectra and light curves of this model which has about 
the same shell mass as Model 3 of lFink et ail (|2010f ) show 
comparably good agreement to observed SNe la as the 
shell-less toy model. 

Taking into account all these results, we argue that 
these systems might yet be promising candidates for 
SN la progenitors. Much more work will be needed to 
properly investigate this possibility. Besides a more de- 
tailed description of the excitation/ionization state in the 
radiative transfer modelling which includes non-thermal 
effects, we need a better understanding of the initial com- 
position of the helium shell and the incomplete burning 
processes which take place in this material to reach re- 
liable predictions of the burning products of the helium 
shell. Although only a tiny fraction of the mass, the 
post-burning composition of the shell material is critical 
to assessing the viability of the sub-Chandrasekhar-mass 
double detonation scenario. 
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